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We study, both analytically and numerically, disorder-induced localization of light in random 
layered structures with magnetooptical materials. The Anderson localization in such structures 
demonstrates nonreciprocal features in the averaged localization length and individual transmission 
resonances. We employ short-wavelength approximation where the localization effects are strong, 
and consider both the Faraday and Voigt magnetooptical geometries. In the Faraday geometry, 
the transmission is strongly nonreciprocal for the circularly polarized waves, whereas in the Voigt 
geometry, the nonreciprocity is much weaker, and it may appear only for the individual transmission 
resonances of the TM-polarized waves. 
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I. INTRODUCTION 

Anderson localization is a well known phenomenon 
associated with wave interference arising from multiple 
scattering by defects. Being originally suggested for the 
suppression of classical diffusive motion of electrons due 
to disorder introduced into a periodic structure [![, the 
Anderson localization is currently associated with many 
physical systems of different origin and nature. As a mat- 
ter of fact, the Anderson localization is universal to all 
wave systems in disordered potentials and media (see, 
e.g., Refs. to cite a few). In particular, the An- 

derson localization occurs upon propagation and scatter- 
ing of light in random media, which are microscopically 
transparent but appear opaque because of the localiza- 
tion effect 0, H, H, 0] ■ Many experimental studies of the 
localization properties of light include the demonstration 
of exponential decay of transmittance with the sample 
length [|| and the study of transverse localization in one- 
and two-dimensional disordered photonic lattices [9|. 

Here we study, both analytically and numerically, the 
Anderson localization of light propagating through ran- 
dom magnetoactive layered structures. We demonstrate 
that an interplay between strong localization and mag- 
netooptical effects produces a number of nonreciprocity 
features in the transmission characteristics. 

Magnetooptical effects and nonreciprocity are wi dely 
exploited in modern optics and applied physics [IcJ 111. 
In particular, magnetoactive periodic structures are cur- 
rently attracting growing attention [l2l . [l3j . The main 
phenomena of interest are the enhanced Faraday effect 
on resonances flij and one-way propagation (nonrecip- 
rocal transmission) (l5l - [l8j employed for the concept of 
optical insulators. The resonant Faraday effect has also 
been shown in connection with the localization of light in 
random layered structures [ll| . Although the destructive 
role of the Faraday effect on the weak localization effect 
in three-dimensional random scattering media was ex- 
amined previously [2(| , there is no analysis of the strong 



Anderson localization in random magnetoactive media. 

In this paper, we examine the transmission proper- 
ties of one-dimensional random layered structures with 
magnetooptical materials. We employ short- wavelength 
approximation, where the localization is strong, and con- 
sider both Faraday and Voigt geometries. In the Faraday 
geometry, magnetooptical correction to the localization 
length results to a significant broadband nonreciprocity 
and polarization selectivity in the typical, exponentially 
small transmission. In the Voigt geometry, averaging 
over random phases suppresses the magnetooptical ef- 
fect, in contrast to the case of periodic structures where 
it can be quite pronounced [3) (3- At the same time, 
in both the geometries we reveal the nonreciprocal fre- 
quency shifts of narrow transmission resonances, corre- 
sponding to the excited localized states inside the struc- 
ture (2ll |22| . This offers efficient unidirectional propaga- 
tion at the given resonant frequency. 

The paper is organized as follows. In Section II we in- 
troduce general formalism for the analysis of the wave 
propagation and localization in random layered struc- 
tures, and also discuss generic aspects of nonreciprocity. 
Section III is devoted to the explicit calculations of the 
averaged localization lengths in the Faraday and Voigt 
magnetooptical structures. Transmission resonances in 
these structures are discussed in Section IV. Finally, Sec- 
tion V concludes the paper. 

II. GENERAL FORMALISM 

A. Basic equations 

We consider transmission of a polarized electromag- 
netic wave incident on a random stack of the length L 
which consists of Af dissipationless magnetoactive layers 
possessing random widths and different optical parame- 
ters. The widths of the layers are independent random 
values with mean value us = L/Af and variance d. Fig- 



ure 1 displays a scheme of the system with two alter- 
nating types of the layers. Here the stack is formed by 
a sequence of Af = 2N layers of two different types la- 
beled by indices 'a' and '6' and the whole structure is 
surrounded by vacuum labeled by index '0'. The z-axis 
is directed across the layers. 

In the simplest case of one type of propagating waves, 
the waves at each point can be described by two am- 
plitudes h v corresponding to the propagation in positive 
v = 1 and negative v = — 1 directions with respect to 
z-axis. Transmission of these waves through the struc- 
ture is described by the total transfer matrix T, which 
expresses the input amplitudes via their output value: 



h(0) = fh(L) 



h = 



(1) 



The amplitudes are assumed to be normalized so that 
the intensity \h\ 2 gives the wave energy. Then the trans- 
fer matrix is unimodular, detT = 1, which ensures the 
energy flux conservation. 

The transmission coefficient T for the wave h + incident 
on the system from the left is simply related to the first 
diagonal element of the transfer matrix: 



T 



1 



T 



(2) 



The corresponding transmittance of the structure is 

t=\t\ 2 . (3) 

Due to Anderson localization, the absolute value of the 
transmission coefficient exponentially decreases with the 
stack length for a typical random realization. Such a 
decrease is described by the dimensionless transmission 
decrement 0,0] 



(InT) 
2N 



(4) 



where (...) stands for the ensemble averaging. In the case 
of Nk 3> I, it coincides with the inverse dimensionless 
localization length £\ oc : 
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The real dimensional localization length is obtained by 
multiplication of £\ oc by the averaged length of one pe- 
riod, 2w. 

We are interested in two transmission characteristics: 
transmittance on a realization, T, Eq. ([3]), and trans- 
mission decrement k, Eq. Following Ref. [7|], we 
represent our system as a sequence of uniform layers of 
media 'a' or '6' and interfaces between them (Fig. [TJ. 
Wave propagation inside the layers is described by diag- 
onal transfer matrices: 



S k = 







(6) 



^1 




FIG. 1: (Color online.) Schematic picture of the wave trans- 
mission and reflection from a random-layered structure con- 
sisting of two types of alternating layers 'a' (here - a mag- 
netoactive material) and '6' (here - air) with random widths. 
Magnetization of the medium, wave polarizations and direc- 
tions of propagation are shown for the Faraday and Voigt 
geometries. 



where ifk is the phase accumulated upon the wave propa- 
gating from left to right through the fc-th layer. Here the 
layers with odd numbers k = 2j — 1, j = 1, 2, N are 
filled with medium 'a', whereas those with even numbers 
k = 2j, j = 1, 2, N are filled with the medium '6'. 

The interfaces are described by unimodular transfer 
matrices F 0a , F ab , F ba , F a0 corresponding, respectively, 
to transitions from vacuum to the medium 'a', from the 
medium 'a' to the medium from the medium '6' to the 
medium 'a', and from the medium 'a' to vacuum. Thus, 
the total transfer matrix (flj of the structure is 



rp E^Oa Q J7 Q Q J7 Q T? Z?& u 
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F 2J _! = F ab , 



F- 



2j 
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F , j — 1,2,..., N. (7) 



Here, in the two last multipliers, we used the group prop- 
erty of the interface transfer matrices: F b0 = p ba p a0 _ 

Our numeric calculations of transmittance ([3|) and 
transmission decrement ^ are based on the exact ex- 
pression ([7]). At the same time, analytical expression for 
the transmission decrement can be readily obtained in 
the short- wavelength approximations. 

First, in the localized regime, the main contribution to 
the localization length is provided by 2N transfer matri- 
ces of the layers and we can neglect the external interface 
transfer matrices F 0a , F a0 just replacing the exact ma- 
trix T by the truncated matrix T': 



T — S1F1S2F2 ... SW-i-FW-iSWF 1 ; 



2JV- 



(8) 
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Second, if the wavelength within the fc-th layer is much 
shorter than the variance of the layer thikness Q: 



< d , 



(9) 



then the phases tpk modulo 2n in the propagation ma- 
trices ^ are nearly uniformly distributed in the range 
(0, 27r). In this approximation, the transmittance for the 
transfer matrix (jSJ) averaged over the phases (fk is re- 
duced to the product of the transmittances of separate 
layers [3| and, furthermore, to the product of transmit- 
tances of the interfaces only Q: 



2iV 



lnCD^lnr, , Tj = l/|(^)n 
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Substitution of Eq. (JTUJ) into Eq. ([5]) yields the follow- 
ing simple expression for the transmission decrement and 
corresponding localization length: 
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This is the main result for calculation of the localiza- 
tion length in the short- wavelength approximation. It 
can be easily extended to any number of alternating lay- 
ers. For instance, considering a random structure con- 
sisting of three types of alternating layers, 'a', and 
'c', with JV = 3N, one has 
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B. Nonreciprocal transmission 

Transmission through a one-dimensional dissipation- 
less linear medium is always reciprocal if there is only 
one (but propagating in two directions) mode in the sys- 
tem. Indeed, while the forward transmission of the h + 
wave incident from the left on the medium is described 
by the 2x2 transfer matrix T, the backward transmis- 
sion of the reciprocal h~ wave incident from the right 
is characterized by the inverse transfer matrix T . Us- 
ing the SU(1, 1) group properties of the transfer matrix, 
one can easily show that the corresponding transmission 
coeffisient ([2]) of the backward wave is T, i.e., exactly 
coincides with that of the forward wave [3, 0| • 

If the system possesses two or more uncoupled modes 
labeled by index a, the situation can be more compli- 
cated. In this case, the waves are marked by the prop- 
agation direction and mode indices: h v ' a . Still, the for- 
ward and backward propagation of each mode a through 
the system (h +,a and hr-" incident waves) are described 
by the 2x2 transfer matrices T a and (T <T )~ 1 character- 
ized by the same transmittance T° . However, the wave 
reciprocal to h + ' a is determined by the time-reversal op- 
eration which changes (because of the k t— > — k 



transformation) but can also affect a [Tl|. In particu- 
lar, if the time reversal operation changes the sign of the 
mode index: a i— > —a, then the reciprocal wave will be 
h' a rather than the backward wave of the same mode, 
/i~ ,<T . Accordingly, the transmittances of the mutually 
reciprocal waves through the system, T* 7 and T~ a , can 
be different. This signals nonreciprocity in the system. 

Note that noreciprocity in the system under consider- 
ation originates from the difference between the modes a 
and — er, and does not depend explicitly on the direction 
of incidence v. Therefore, in practice, it is sufficient to 
compare only forward transmissions of the modes ±er, de- 
scribed by the transfer matrices T ±cr and transmittances 

As we show below, propagation of light in magne- 
toactive layered media offers nonreciprocal transmission 
which can be explained within the above formalism. 
There are two main geometries typical for magneto- 
optical problems (lo| : (i) the Faraday geometry, where 
the magnetization is collinear with the direction of propa- 
gation of the wave, and (ii) the Voigt (or Cotton- Mouton) 
geometry, where the magnetization is orthogonal to the 
direction of propagation of the wave (see Fig. 1). In the 
next Sections we study the averaged transmission decre- 
ment (Sec. Ill) and individual transmission resonances 
(Sec. IV) in both geometries. 



III. LOCALIZATION DECREMENTS 

A. Faraday geometry 

In the Faraday geometry both magnetization and the 
wave vector are directed across the layers, i.e., along the 
z-axis (see Fig. 1). We assume that the magneto-optical 
effects are described exclusively by the dielectric tensor, 
while the magnetic tensor is equal to one. In this case, 
the electric induction in the medium reads lCf 



D = eE + ig x E = eE 



(0,0, Q), 



(13) 



where E is the electric field of the wave, e is the isotropic 
dielectric constant in the absence of magnetization, and 
Q = eq is the magneto-optical constant proportional to 
the magnetization of the medium. In what follows, we 
assume that \q\ <C 1 and will be interested in the effects 
linear in q. 

Thus, the dielectric tensor in the Faraday geometry 
has the form: 



(14) 



Solving stationary Maxwell equations for the magnetic 
field H = {H x , H y , 0) in a homogeneous medium, 




k x [e- 1 (k x H)] 



fcoH 



(15) 
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fco = w/c, we find that the eigenmodes of the problem 
are circularly polarized waves: 



V2 




v,a = ±1. 



The corresponding wave electric field E is 



E' 



'" u TTV .(7 

iva—il ' . 

k 



(16) 



(17) 



In these equations, H v,ry are the wave amplitudes, 
whereas 



k = nfco\/l + 



aq, 



(18) 



is the propagation constant affected by the magnetization 
parameter q and depending on a. In the linear approxi- 
mation in q, k ~ nko(l + aq/2). 

Parameter a is the mode index which determines the 
direction of rotation of the wave field. In this manner, 
the product vo~ represent the hclicity 



(19) 



which distinguishes the right-handed (x = +1) and left- 
handed (\ — —1) circular polarizations defined with 
respect to the direction of propagation of the wave. 
Note that the time reversal operation keeps helicity un- 
changed, whereas a changes its sign 11]. Thus, the re- 
ciprocal wave is given by H~' J, ~ <T , precisely as described 
in Section IIB. 

Consider the wave transformation at the interface be- 
tween the media 'a' and The helicity of the wave flips 
upon the reflection and remains unchanged upon trans- 
mission. As a result, parameter a remains unchanged, so 
that there is no coupling between the modes with a = +1 
and <7 = —1 (see Fig. 1), and these modes can be studied 
independently. From now on, for the sake of simplicity, 
we omit a in superscripts and write explicitly the values 
of the direction parameter v = ±1. In this manner, the 
boundary conditions for the wave electric and magnetic 
fields at the 'a'-'b' interface read 
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Substituting Eqs. and (HTJ) into Eqs. $2U\), we ob- 
tain that the wave amplitudes in the two media, H a b = 

\ T 

H^ b , H~A , are related through the unnormalized in- 
terface transfer matrix: 



H a = F ab H b 



1 



F ab — 



kb + k a kb k a 



2k h \ h-ka k b + k a 



(21) 



Here k a ,b are the wave numbers (|18p in the corresponding 
media. 

The determinant of the matrix (f2"Tj). deiF ah — k a /k 0: 
determines the choice of the normalized amplitudes 
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FIG. 2: (Color online.) Localization decrement k vs. 
magneto-optical parameter Q for opposite modes propagating 
through a two-component random structure in the Faraday 
geometry (see details in the text). The modes with a = ±1 
correspond to either opposite circular polarizations or prop- 
agation directions. Numerical simulations of exact equations 
(symbols) and theoretical formula (|24|) (lines). 



and the normalized interface transfer matrix is 
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kb k a 
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(23) 



2^fk a k b 

Considering now a random multi-layer structure and 
calculating the localization decrement from Eq. (TTT|) with 
Eqs. (|18|) and (|23l) . we obtain in the linear approximation 
in q: 



2 In- 



„(o )=ln K + ^) 2 

An a nb 
(1) cr ,n a 

2 n a + n b 



rib 



(24) 



Thus, the localization decrement acquires the first-order 
magneto-optical correction caused by the Faraday 
effect. This correction depends on tr, i.e., on the polar- 
ization helicity x an d the- propagation direction v through 
a = x v - F° r the- reciprocal waves with the same x an< i 
opposite v, has opposite signs. This signals nonre- 
ciprocal localization in a Faraday random medium. In 
practice, the nonreciprocal difference in the transmission 
decrements (|2~4"|) can be observed by changing sign of ei- 
ther propagation direction v (with the helicity fixed) , or 
polarization x, or magnetization q. 

Despite the magneto-optical correction to the localiza- 
tion decrement is small in magnitude, k^°\ it 
still might result in a significant difference in the typi- 
cal transmission spectrum. This difference is described 
by an additional factor of oc exp[— 2Nk^] in transmit- 
tance, which is exponential with respect to the length 
of the structure. Hence, small correction (|24"|) brings 
about significant broadband nonreciprocity or polariza- 
tion selectivity in the typical small transmission when 
N\k^\ > 1. 



5 



Fig. [2] shows dependence of the localization decrement 
on the magnetization parameter Q = eq calculated nu- 
merically and compared to analytical result (|24|) . Numer- 
ical simulations were performed for the structure contain- 
ing Af = 2N = 90 alternating layers of air (e = 1, Q = 0), 
and bismuth iron garnet (BIG), with dielectric constant 
e = 6.25 and magneto-optic parameter reaching Q — 
0.06. The thicknesses of layers were randomly distributed 
in the range 50-^150 /im (i.e., w = 100/^m, d = 50/mi), 
whereas the excitation wavelength was 632 nm. The aver- 
aging was performed over 10 5 realizations of the random 
sample. One can see excellent agreement between nu- 
merical simulations and analytical results showing linear 
splitting of the <j = 1 and <j = — 1 localization decrements 
as a function of the magneto-optic parameter. 



B. Voigt geometry 

In Voigt geometry the magnetization is directed or- 
thogonally to the z-axis, say, along the y-axis. Then, the 
dielectric tensor is 




(25) 



The first-order interaction of the wave with the magneti- 
zation occurs only upon oblique propagation of the wave 
in the xz-plane, i.e., when k x — const ^ (see Fig. 1). 
The eigenmodes in such problems are the TE and TM 
linearly-polarized modes. Since the TE mode is uncou- 
pled from the magnetization and effectively propagates 
as in isotropic layered medium, we consider only non- 
trivial TM polarization. Solving Maxwell equations (fTS"|) 
with dielectric tensor (1231) we obtain that the wave mag- 
netic field of the TM mode is directed along the y-axis 
whereas the electric field has x- and z-componcnts. For 
transmission and localization properties of the structure, 
only components tangential to the layer interfaces are 
important: 
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Here parameters v = ±1 and er = ±1 indicate propa- 
gation in the positive and negative z and x directions, 
respectively, k\\ — y/k 2 — t2 1 — I wliorooo 
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, k = nWl - Q 2 - (27) 



Note that A v ' a ~ (ufeii + iaqk±) / (eko) and k ~ nko in 
the linear approximation in q, so that the magnetization 
affects imaginary parts (i.e., phases) of the amplitudes 
A v,a and does not affect the propagation constant, cf. 
Eqs. (|T7| and ffT5|). 

In the Voigt geometry, direction of the transverse wave 
vector component, er, serves as the mode index. The 
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FIG. 3: (Color online.) Localization decrement k vs. 
magneto-optical parameter Q for reciprocal waves with a 
: I propagating through a three-component magnetooptical 
random structure in the Voigt geometry (see details in the 
text). Numerical simulations of exact equations (symbols) 
and theoretical formula (|32|) (line). 



mutually reciprocal waves are H v,a and H~ V >~ C! because 
the time reversal transformation reverts the whole wave 
vector, k H> — k. Thus, we again deal with a formalism 
described in Section IIB. 

Evidently, the parameter a is not changed upon reflec- 
tion and transmission through the layers, i.e., the modes 
with a = ±1 are uncoupled from each other. Therefore, 
for the sake of simplicity, we omit the mode index in 
superscripts, and write explicitly only the values of the 
direction parameter v — ±1. Matching the tangential 
components of the fields ([26)) at the interface between 
media 'a' and Eq. (|20j) . we obtain the unnormalized 
transfer matrix relating the wave amplitudes in two me- 
dia, H„ = F ab H h M: 



F 



I 



A^ A„ 



A b A„ 



(28) 



Here A^ b are the amplitudes (|2~T1) in the corresponding 
medium. Noticing that A~ = A + * , we calculate the de- 
terminant of the matrix <j28j), det^ afc = ReA£/ReA+, 
which determines the normalized field 



V2ReI+ H a , h b = J 2RcA+ Hi 



lf) ,i 6 . (29) 

As a result the normalized transfer matrix takes the form 
1 



F = 



AReA+ReAi 



A+ I A+* A+* _ A+-- 

At-A+ A+ + A+* 



.(30) 



In contrast to the Faraday geometry, in the Voigt 
geometry the linear magneto-optical correction changes 
only phases of the transmission and reflection coeffi- 
cients, whereas corrections to the interface transmittance 
start with the terms oc q 2 . In short-wave limit, only 
these transmittances determine the total transmittance, 
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FIG. 4: (Color online.) Transmission spectra of a random 
magneto-optical sample in the Faraday geometry (see details 
in the text) for waves with a = ±1. While the averaged 
localization decrements are only slightly different (Fig. 2), all 
individual resonances are shifted significantly as compared 
with their widths, Eq. (|33[) . 



Eq. (fl"U|). Therefore, a short- wavelength transmission 
through a random multilayered stack is reciprocal and 
is not affected by magnetization in the first-order ap- 
proximation. In particular, substituting Eqs. (1271) and 
(1501) into Eq. (ITT1) . we arrive at the localization decre- 
ment for a two-component random layered structure (cf. 
Ref. (H): 



k = In 



(eifc 2 || + e 2 fci|| 
4eie 2 fci||fc2|| 



+ 0(q 2 ), 



(31) 



where ~ y/ n 2 fcp — A; 2 . Obviously, it depends on nei- 
ther propagation nor magnetization directions. 

It is worth remarking that transmission through a pe- 
riodic structure consisting of the two types of alternat- 
ing layers is also reciprocal in the Voigt geometry for 
all wavelengths. However, a periodic structure with a 
cell consisting of three different layers (which breaks the 
mirror reflection symmetry) can demonstrate significant 
nonreciprocity [lg, [l8j]. At the same time, the short- 
wavelength localization in a random-layered structure is 
still reciprocal for any number of components. This is 
because of independent action of all the interfaces after 
the phase averaging, Eq. (jlOp . In particular, considering 
a random structure consisting of three types of alternat- 
ing layers, 'a', '6', and 'c', we substitute Eqs. (|27|) and 
([50)1 into Eq. (TT2"j) and obtain the localization decrement 
in the first-order approximation in q: 



hi 



(gqfc||6 + £fefc||q) (gbfc||e + fcfe||b) (gcfc||a + Egk\\c) 

8£ a e 6 e c fc|| Q fc|| b fc|| c 

(32) 

We verified this result numerically as it is shown in 
Fig. [3] There we calculated forward and backward, v = 
a = ±1, transmission through M — 3N — 60 random- 
width layers of three alternating types: air, BIG, and 
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FIG. 5: (Color online.) Differential transmittance, T + — 1 , 
for two resonances from Fig. 4 as dependent on the value of 
magneto-optical parameter Q, cf. Eq. (1331) . 



glass (n = 1.5). The excitation wavelength was 632 nm 
and the angle of incidence #o = sin (fc^/feci) = 23°. The 
layer thicknesses were randomly distributed in the range 
50-^150 /im (w = 100/im, d = 50/mi), and the averaging 
was performed over 105 realizations of the sample. The 
results of simulations indicate no magnetooptical effect 
and completely reciprocal localization. 



IV. TRANSMISSION RESONANCES 

Averaged localization decrement is associated with ex- 
ponential decay of the incident wave deep into the infi- 
nite sample 01 • For a finite sample, this is so only 
for typical realizations. However, there exist some reso- 
nant realizations of the sample at a given frequency (or, 
equivalently, resonant frequencies for a given realization) 
where transmission is anomalously high and is accompa- 
nied by the accumulation of energy inside the sample. 
[2ll . [22| Such resonant transmission corresponds to ex- 
citation of the Anderson localized states (quasi-modes) 
inside the sample. Akin to the resonant localized states 
in photonic crystal cavities, the transmission resonances 
in random structures are extremely sensitive to small per- 
turbations: absorption, 22| nonlinearity, [24| and, as we 
show here, magnetoactivity. 

Figure |4] shows transmission spectra for two modes 
(7 = ±1 (i.e., either with opposite helicities or propa- 
gation directions) in one realization of a magnetooptical 
sample in the Faraday geometry. The parameters of the 
sample are the same as in Section IIIA with Q = 0.06. 
One can see strong splitting of the a = ±1 transmis- 
sion resonances which have exponentially narrow widths 
oc Kexp(—K,N)/2w [22]. This offers strongly nonrecip- 
rocal, practically unidirectional, propagation or polariza- 
tion selectivity in the vicinity of resonant frequencies. 

To estimate the splitting of resonances, we note that 
the wavenumbers in magnetooptical materials are shifted 
due to the Faraday effect, Eq. (fP5|) . Hence, the shifts 
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FIG. 6: (Color online.) Differential transmittance, T + — 7 , 
for in the vicinity of a single resonance in the Voigt geometry 
(see Section IIIB for details) as dependent on the magneto- 
optical parameter Q. 



of the resonant wavenumbers of the random Faraday 
medium can be estimated by averaging of this shift over 
different materials in the structure: 



Afc r , 



qnk 



(33) 



where (...) stands for some average of (...). Using qn ~ 
{Qa n a + qb n b)/2 f° r estimation in the two-component 
structure, we obtain AA rcs ~ — a 3.6 nm, which agrees 
with the cr-dependent splitting observed in Fig. [4] 

Figure [5] displays the differential transmission for the 
waves with a = +1 and a = —1 as a function of magneto- 
optical parameter Q for two resonances lying in a narrow 
frequency range in Fig. [4] In agreement with estimation 
(|33| , one observes the linear dependence of the resonance 
splitting on magnetization. 

In the Voigt geometry, the resonances also allow non- 
reciprocal transmission and demonstrate splitting of the 
resonant frequencies. In Fig. [6] we show differential trans- 
mission for reciprocal waves with a = ±1 in the vicinity 
of one resonance for the three-component structure con- 
sidered in Section IIIB. The splitting is very small in 
this case, and a = +1 and a = — 1 resonances overlap 
significantly. Because of this, the differential transmit- 
tance in Fig. 6 is tiny, its amplitude linearly grows with 
Q, whereas the frequency positions of its maximum and 
minimum correspond to the width of the original reso- 
nance and are practically unchanged. 

Unlike the wave-number shift in the Faraday geome- 
try, the noreciprocal shift of resonant frequencies in the 
Voight geometry arises from the phases of the amplitudes 
A, Eq. (|27|) . These phases are responsible for the phases 
of transmission coefficients between the layers and can 
be estimated as ~ q(akx)/(vk») = qta,n9, where 9 
is the angle of propagation with respect to the z-axis. 
The phases accumulated at a layer effectively shift the 
wave numbers as uAfcii = Afccos6* ~ S/w, where w is 



the thickness of the layer. Averaging over different ma- 
terials in the random layered structure, we estimate the 
nonreciprocal shift of the resonant wave number: 



Afc r , 



q sin o 
w cos 2 9 



q\ sin#| 
w cos 2 9 



(34) 



This shift is cr-dependent, i.e., nonreciprocal, and much 
smaller than the Faraday-geometry shift (|3"3")l as kw > 
kd ^> 2-7T in the short-wavelength limit, Eq. ©. For 
the parameters in use, with Q = 0.06, we have AA rcs ~ 
— cr3 • 10~ 4 nm, which agrees with the data plotted in 
Fig. 6. 



V. CONCLUSIONS 

We have studied the transmission and localization of 
light in magnetoactive layered structures. An interplay 
between the Anderson localization and magnetooptical 
effects brings about various nonraciprocal phenomena in 
the transmission characteristics. We have analyzed the 
effects of the medium magnetization on the wave trans- 
mission in both Faraday and Voigt geometries in the 
short- wavelength limit. 

Specifically, in the Faraday geometry the averaged lo- 
calization length acquires the first-order magnetooptical 
corrections of the opposite signs for the opposite prop- 
agation directions or opposite circular polarizations of 
light. This leads to a broadband nonreciprocity or po- 
larization selectivity in the typical exponentially small 
transmission observed in such structures in the regime 
of the Anderson localization. At the same time, random 
transmission resonances acquire significant nonreciprocal 
frequency shifts which result in efficient unidirectional 
propagation at the given resonant frequency. In the Voigt 
geometry, for the TM-polarized waves, the localization 
length is always reciprocal in the first-order approxima- 
tion, whereas the transmission resonances show nonre- 
ciprocal frequency shifts but much smaller than those in 
the Faraday geometry. 

Thus, we have observed that disorder-induced localiza- 
tion of light in random layered structures with magne- 
tooptical materials demonstrates nonreciprocal features 
in both the averaged localization length and individual 
transmission resonances. Our results demonstrate that 
the Anderson localization can significantly enhance the 
magnetooptical effects, and this property can be em- 
ployed for a design of novel types of efficient nonrecipro- 
cal devices which do no require periodicity and specially 
designed cavities. 
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